Dissipation and criticality in the lowest Landau level of graphene 
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The lowest Landau level of graphene is studied numerically by considering a tight-binding Hamil- 
tonian with disorder. The Hall conductance a^y and the longitudinal conductance o-:,^:x. are computed. 
We demonstrate that bond disorder can produce a plateau-like feature centered at i^ = 0, while the 
longitudinal conductance is nonzero in the same region, reflecting a band of extended states between 
±Ec, whose magnitude depends on the disorder strength. The critical exponent corresponding to 
the localization length at the edges of this band is found to be 2.47±0.04. When both bond disorder 
and a finite mass term exist the localization length exponent varies continuously between ~ 1.0 and 
~7/3. 



The quantum Hall effect (QHE), either integer or 
fractional, is one of the most intriguing phenomena in 
physics [T] that has drawn, and still continues to draw, 
enormous attention even after two decades since its dis- 
covery. In graphene unconventional QHE [21 |3] corre- 
sponding to CTxy = ±(n-|-l/2)4(e^//i) = ±4i/(e^//i), n = 
0, 1, 2, . . ., where a^y is the Hall conductance, has been 
addressed in terms of low lying excitations akin to rela- 
tivistic Dirac fermions [H [5] ; the factor of 4 arises from 
the two- fold valley and spin degeneracies; e is the elec- 
tronic charge and h the Planck's constant. Yet, a more 
complete theory of QHE requires an understanding of the 
localization-delocalization transitions within the Landau 
levels, which reflects a very special quantum phase tran- 
sition, especially in graphene with its low energy Dirac 
spectra O [7l |8]. It is indeed paradoxical that such a 
precise quantization requires material defects and disor- 
der, and even the nature of disorder seems to matter for 
graphene. 

A remarkable recent discovery in graphene is, what 
appears to be, a z^ = plateau in a sufficiently large 
magnetic field [SJ [TU] , although the plateau in the mea- 
sured pxy is difficult to decipher. The longitudinal resis- 
tivity pxx, on the other hand, is observed to be nonzero 
(> /i/e^) in the region of the plateau |TD], in sharp con- 
trast to the non-dissipative behavior of pxx in conven- 
tional QHE, crying out for a theoretical explanation. It 
is most peculiar because "a dissipative quantum Hall 
plateau" is an oxymoron, for, according to Laughlin fl] , 
the precise quantization requires zero dissipation. An 
intriguing explanation of this paradoxical phenomenon 
is given in Ref. |10j . where the removal of spin degen- 
eracy plays a special role, and the nonzero longitudinal 
resistivity pxx is ascribed to a pair of gapless counter- 
propagating chiral edge modes carrying opposite spins, 
while a spin gap in the bulk protects the quantization of 
the Hall conductance. The role of interactions is crucial 
in this theory. 

Surprisingly, in the present Letter we shall demon- 
strate by explicitly computing a^y, and the longitudinal 
conductance axx , that the existence of a plateau- like fea- 
ture at z^ = and the non-zero p^x can be attributed 



to a single mechanism, an inter-valley coupling that is 
induced by bond disorder. Because Gxy is zero at u — Q, 
Pxx = l/fxx at the same point. The most remarkable fea- 
ture is that there is a band of extended states centered 
at zero energy, whose extent depends on the strength of 
bond disorder. In contrast to Ref. [10], in our picture 
there is no quantization of axy at :^ = 0, only a sloping 
behavior as a function of energy. In addition, dissipation 
is not an edge phenomenon, but a bulk one. 

Another remarkable aspect of the integer QHE in 
graphene is the possible existence of a continuously vary- 
ing critical exponent of the divergence of the localization 
length within the Landau band for a class of disorder. We 
show: (1) If there is only bond disorder, the critical expo- 
nent is close to ~ 7/3, same as in the conventional integer 
QHE; (2) when a finite mass term is added in addition 
to bond disorder, the critical exponent depends on the 
ratio of the intensity of bond disorder to the finite mass, 
continuously varying from ^ 1.0 to ^ 7/3, as the system 
is tuned from weak to strong disorder, which correctly 
reaffirms the results of Ref. ^ obtained by an entirely 
different method; the present method is more powerful, 
because larger system sizes can be handled. A finite mass 
appears to be experimentally relevant. Although the re- 
moval of the nodal degeneracy in the lowest Landau level 
may be of many body origin, it can be approximately ac- 
counted for by including an explicit mass term in the 
Hamiltonian. 

We study the tight-binding model of graphene sub- 
ject to a constant perpendicular magnetic field and dis- 
order [111 , using real space transfer matrix and exact 
diagonalization methods, to explore the peculiar prop- 
erties of the lowest Landau level discussed above. The 
Hamiltonian defined on a honeycomb lattice of dimen- 
sion Lx X Ly, shown in Fig. IT] is 



H ==^[(e„,A + ■m)cl^j^Cn.A + (e„,s - ''n)clgCn,B] 

n 

3 

- Yl Yitn,ke'''"''cl^Cn+S„B + h.c), 
n k = l 



(1) 



where the summation of n ranges over all unit cells, and 
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FIG. 1: (Color online) Graphene lattice. Lattice sites belong- 
ing to the same shaded horizontal or vertical stripes share 
the same indices i or j. The solid circles correspond to the 
A sublattice and the open circles to the B sublattice. The 
three nearest neighbor vectors joining the two sublattices are 
6k, k = 1...3. 



Cn,Ai Cn,_B are the fcrmionic annihilation operators in the 
unit cell n for the sublattices A and B, respectively. The 
spacing between vertical slices is v3a/2, where a is the 
bond length. In the transfer matrix calculations a is the 
unit of length. The size of the sample is chosen such 
that L^ = N{V3a/2) and Ly = M(3a), where TV is the 
number of vertical slices and M is the number of unit 
cells on a vertical slice. 

The spin degrees of freedom are omitted, as we assume 
that the magnetic field is sufficiently strong to completely 
polarize them. The on-site energies £n.A and en,B are in- 
dependent random variables. Thus, Vn = (en,A + en,B)/2 
is a random potential and M„ — (en, a ~ £n,B)/2 is mass 
disorder in the corresponding language of the low energy 
spectra of Dirac fermions. Here we do not consider mass 
disorder and therefore choose Cn.A = Cn.s = Vn, where 
Vn is uniformly distributed in the range [—gv/'2,gv/'2]- 
The mass m provides a charge density modulation be- 
tween the two sublattices and leads to an energy gap. In 
the Dirac language this gap would appear as a parity- 
preserving mass. We choose tn.k = t + Atn,k and set 
t — I, providing a natural energy scale. The quan- 
tity Atn,k is a random variable uniformly distributed 
in the range [—gx/2,gx/2\, characterizing the bond dis- 
order. The phases an,k are such that the magnetic 
flux per hexagonal plaquette, (j), is 1/Q, in units of the 
flux quanta 4>o = h/e. We choose a gauge such that 
On,! = T^i/Q for the vertical bonds in slice i as in Fig. IT] 
and an_2 = On.s = 0. 

The longitudinal conductance a^x is studied using 
the well developed transfer matrix method. Consider 
a quasi- ID system, L^ 3> Ly with a periodic bound- 
ary condition only along the y direction. Let ^^ = 
{4'i,i,'4'i,2, ■ ■ ■ ,'4'i,2M)'^ be the amplitudes on the slice i 
for an eigenstate with a given energy E; then amplitudes 
on the successive slices are related by the matrix multi- 



plication: 






1 



*, 
^^-1 



, (2) 



where 7^ is a diagonal matrix with elements 
{ti^i,ti^2, ■ ■ ■ ,ti.2M) representing the hopping matrix 
elements connecting the slices i and i + I, and Hi is the 
Hamiltonian within the slice. All postive Lyapunov ex- 
ponents of the transfer matrix [T3], 71 > 72 > • ■ • > 72 Af, 
are computed by iterating Eq. (p| and performing 
frequent orthonormalizations. The convergence of this 
algorithm is guaranteed by the well known Osledec 
theorem [TJ. The conductance per square, a^x, is given 
by the Landauer formula [T^ [T51 [T71 [T5] (note the special 
factor of \/3 in the argument of cosh) : 



2M 

E 



1 



T^cosh2(2V3M7,)' 



(3) 



The localization length in the quasi-lD system of width 
Ly is given by Am = 1/72M- Assuming single parameter 
scahng, Xm/M = f{\E - Ec\M^''''), the data collapse 
yields the critical exponent vi and the critical energy E^. 
To compute the Hall conductance a^y, we impose pe- 
riodic boundary conditions in both directions of the sys- 
tem. The Hamiltonian (IT|) is diagonalized to obtain a 
set of energy eigenvalues E^ and the corresponding set 
of eigenstates |a) for a — 1, . . . , 2M x N 
computed using the Kubo formula |12j : 



Then a^y is 



^xy 



{E) 



LxLy 



E 



Eo,<E<Eg 



(a\vy\f3){f3\vx\a) - {x ^ y) 
(E^ - EpY 



(4) 
where Vx = [H,x]/ih is the velocity operator along the 
X direction and similarly for Vy in the y direction. Note 
that the bonds in Fig. [T] that are not parallel to the y 
direction contribute to both v^ and Vy. The summation 
corresponds to sum over the states below and above the 
energy E. Finally, the expression is disorder averaged. 

In the language of Dirac fermions [7] random hopping 
gives rise to both intranode and internode scattering be- 
tween states on different sublattices. Intranode scatter- 
ing appears as a random abelian gauge field, and the two 
inequivalent nodes have opposite charges corresponding 
to this gauge field. When projected to the lowest Lan- 
dau level, the abelian gauge field leaves it unaffected. 
However, the internode scattering mixes the degenerate 
states corresponding to the two inequivalent nodes and 
produces extended states at zLEc- The existence of ex- 
tended states at energies symmetric about E = is the 
consequence of the sublattice symmetry of the disorder 
(often referred to as the chiral or the particle-hole sym- 
metry). It is this special symmetry that leads to a di- 
vergent density of states and delocalized states a,t E = 
[71 [IH]. However, the calculated finite axx and a linear 
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FIG. 2: (Color online) The longitudinal conductance (Txx as a 
function of energy E. The magnetic flux through the hexag- 
onal plaquette (fi = 1/200 and qt ~ 0.5. 
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FIG. 3: (Color online) Longitudinal conductance CTxx a.t E = 
as a function of the transverse size Af for gr ~ 0.5 and 
gv ~ 0. The inset shows a^^{E = 0) as a function of qt for 
M = 32. The value appears to be independent of g-r and 

2 

close to -^ ; the errors bars are about the size of the scatter 
in the numerical data. The near universality disappears if gv 
is added. 



variation of a,j.y with a small slope in the energy range 
~Ec < E < Ec hint at the existence of a band of delo- 
calized states between ±Ec, as shown below. 

In the transfer matrix calculation of CTxx (see Fig. l2l, 
we chose a magnetic field <j) = 1/200. An iteration of 
Eq. (|2| of the order of 10^ to 10^ was performed until 
the relative errors of less than 0.5% of all the Lyapunov 
exponents were achieved. The longitudinal conductance 
CTxx at exactly E ~ is computed according to the Lan- 
dauer formula, Eq. ([3| , in systems with different values of 
M, as shown in Fig. |3] For only bond disorder, gx — 0.5, 
a non-zero longitudinal conductance is observed for all 



dependent of M and gr- This behavior of a^x implies 
an unusual dissipative nature. We have also checked the 
existence of the dissipative behavior when g^ 7^ in ad- 
dition. 



1.0 



(a) 



1.0 



0.1 




■ M=16 

• M=32 

, ▲ M=48 

*V. ▼ ^^^^ 

^^ ^ M=80 



0.01 ,-0.02 



0.03 



0.1 
0.001 



0.010 



0.100 



2 




M"ME-EJ 








(b) . 


■•■-., 




■ M=16 
• M=32 
▲ M=48 


1.0 

> 




■\ 


♦4 


▼ M=64 
^ M=80 

* 



M"'' (E-E 



0.001 



2.5 



2.0 



1.5 



1.0 



(C) 



u,^— * 



0.10 



St 



1.00 



system sizes. The value of cr-x 



is found to be in- 



FIG. 4: (Color online)(a) Scaling curve for the case gr = 0.5 
and 771 = 0. Insert shows the size dependence of the local- 
ization length computed in systems with different sizes, (b) 
Scaling curve for the case gr ~ 0.2 and m — 0.2. (c)The criti- 
cal exponent i/g as a function of the random hopping intensity 
gr for m = 0.2. It varies continuously from ve = 1.10 ± 0.05 
to 1/ = 2.36 ± 0.07 as the disorder intensity gr is tuned from 
the weak to strong. The dashed line corresponds to ue = 7/3. 



To the study of the critical behavior in the presence 
of bond disorder in the massless case, we set gr = 0.5. 
The critical exponent is expected to be independent of 
the value gr [Z]. The renormalized localization lengths 
Xm/M as a function of E for various M are plotted in 
in the insert of Fig. [IF a). The critical energy Ec is lo- 
cated at a non-zero value Ec — 0.0167 where Xm/M is 
independent of M. A successful data collapse based on 
the data with E > Ec leads to a critical exponent of 
1/ = 2.47 ± 0.04, close to conventional integer QHE. The 



scaling form is depicted in Fig. [IF a). 

For the massive case, we vary bond disorder with a 
fixed m = 0.2 for the purpose of illustration. An example 
of data collapse is shown in Fig. l4|b) with gx = 0.2. 
The critical exponent vi — 2.08 ± 0.01 is different from 
conventional QHE; vi varies continuously from 1.1 ±0.05 
to 2.36 ± 0.07 as the system is tuned from gx = 0.02 
to gT = 1-5. This behavior agrees with the previous 
results [?]■ When g^ ^ m, the effect of the finite mass 
is negligible, and the critical exponent of vg ~ 7/3 is 
recovered, as before. 
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FIG. 5; (Color online) Hall conductance a^y as a function 
of energy E with bond disorder intensities gr- The i/ — 
plateau emerges due to the bond disorders. 

In computing the energy dependence of Uxy a relatively 
large flux <j> = 1/20 is chosen, because smaller values of 
flux involve many Landau bands in the diagonalization 
calculations, which are hard to track accurately. The 
chosen system size was N = M = 40, and an average over 
1000 disorder realizations was performed. The results are 
shown in Fig. [5) Although there is not a strict plateau 
at J^ = 0, there is a break in the slope in the rise between 
ly = —1 toi^= 1, as the energy sweeps past the band 
center at E = 0. This can be construed as a plateau- 
like feature. Since the lowest Landau level splitting is 
expected to increase with gT, so is the extent of the region 
around E = with a smaller slope. 

The striking results here are the band of extended 
states in the region —Ec < E < Ec for bond disorder 
and the vindication of a continuously varying localiza- 
tion length exponent when, in addition, there is a finite 
uniform mass present in the Dirac spectrum. The situ- 
ation is a bit more subtle, however. We had previously 
observed that the density of states diverges very weakly 
at i? = [7 . In particular, for the Lorentzian distri- 
bution of disorder this divergence was found to be ex- 
actly logarithmic. A log E divergence was predicted in 
Ref. 112] for Gaussian disorder corresponding to a Hamil- 
tonian which in fact is formally identical when projected 
to the lowest Landau level. Such a weak divergence at 



E — may give rise to fluctuations responsible for dissi- 
pation leading to a flnite Uxx [3 [H] . The slightly sloping 
profile of (Txy centered at i? = is harder to explain ana- 
lytically. The dissipative behavior at z/ = is consistent 
with experiments. However, in contrast to Ref. [TD], this 
dissipative behavior is a bulk phenomenon, not an edge 
phenomenon. It is possible to test this experimentally 
by varying the aspect ratio of the sample. But we only 
have a sloping plateau-like feature of a^y at ;/ = un- 
like Ref. [T^. It is possible that the difference between 
the two pictures depends on the relative size of the spin 
splitting compared to the width of the extended band of 
states. Clearly further experimental and theoretical work 
would be very helpful to elucidate the precise nature of 
this exciting new development. 
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